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dau two-fluid hydrodynamics, we predict a damped mode centered at zero 
frequency, in addition to the usual second sound doublet. 
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I. INTRODUCTION 



Since the original discovery in 1995, the subject of Bose-Einstein condensation in trapped 
atomic gases has become a major field of research Even though these systems consist of 
a very dilute vapor, they exhibit robust collective oscillations which are strongly infiuenced 
by mean-field interactions and collisions between the atoms. There has been considerable 
theoretical work devoted to describing the collective modes of a trapped gas at very low 
temperatures in terms of the solution of the time-dependent Gross-Pitaevskii (GP) equation 
for the macroscopic wavefunction of the condensate. As discussed in several recent reviews 
|T|J^, there is excellent agreement between experimental observations for T ^ Tbec and 
theoretical calculations based on the T = GP equation. 

At elevated temperatures where the condensate is appreciably depleted by thermal ex- 
citations, one enters a more complex regime in which collisions between the atoms must be 
considered. Two limiting cases for the dynamics of the gas correspond to the collisionless 
and hydrodynamic regimes In the collisionless regime, the main effect of the non- 

condensate component appears to be a shift in the collective mode frequencies as a result 
of the change in the condensate number, and to the appearance of Landau damping. There 
have been several studies of this collisionless region in trapped gases where dynamic mean 
fields dominate the physics (we refer to the review articles in Ref. 0). The second regime 
arises when collisions between atoms are rapid enough to establish a state of dynamic local 
equilibrium in the non-condensate gas. This collision-dominated hydrodynamic regime is 
the subject of this paper. To be in this regime the collective modes of frequency u must 
satisfy the condition ujt -C 1, where r is some appropriate relaxation time for reaching local 
equilibrium. As a rough estimate, we can take r to be the collision time for a classical gas. 
Here l/r^ = V2 nav, where a = Svra^ is the low energy quantum mechanical cross-section 
for bosonic atoms and v is the average thermal velocity |^. It is clear that a high density 
and/or a large atomic scattering cross-section are favorable for reaching the hydrodynamic 
region. Experiments which can probe this region in trapped Bose gases are now feasible and 
promise to provide much new physics. 

One finds in the hydrodynamic regime that the dynamics of the condensate and non- 
condensate components can both be described in terms of a few macroscopic (coarse-grained) 
variables (such as the local densities and velocities of the two components). The coupled 
equations of motion for these local quantities will be referred to as the two-fiuid hydrody- 
namic equations. The microscopic basis of these two-fiuid equations rests on a generalized 
Gross-Pitaevskii equation for the condensate atoms and a quantum kinetic equation for the 
non-condensate atoms. These two components are coupled through mean-field interactions 
as well as collisions between the atoms. The authors have recently given a detailed deriva- 
tion and discussion of such a two-fiuid hydrodynamics for trapped atomic gases at finite 
temperatures 

These equations were derived for temperatures where the dominant thermal excitations 
in the trap can be treated as atoms moving in a self-consistent Hartree-Fock field. In 
particular, the original ZGN hydrodynamic equations derived in Ref. were generalized in 
Refs. 0,^ to include collisions between condensate and non-condensate atoms. This allows 
for the possibility of treating the situation in which atoms of the non-condensate are in local 
thermodynamic equilibrium among themselves but are not in diffusive equilibrium with the 
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condensate atoms. The resulting ZGN' hydro dynamic equations involve a characteristic 
relaxation time which is the time scale on which local diffusive equilibrium is established. 
This equilibration process leads to a novel damping mechanism which is associated with the 
collisional exchange of atoms between the two components. The ZGN' equations are briefly 
reviewed in Section II. 

In Section III, we further generalize the ZGN' equations by considering the effects of 
deviations from local equilibrium within the non-condensate. At the finite temperatures 
of interest, this deviation from local equilibrium gives rise to damping associated with the 
thermal conductivity and the shear viscosity. This generalization has already been discussed 
in Section V of Ref. |Tl|] starting from the ZGN hydrodynamic equations. 



In the limit that the two components are in complete local equilibrium with each other, 
our two-fluid hydrodynamic equations reduce to those first derived by Landau in 1941 [0. 
The Landau two-fluid equations give an excellent description of the low frequency response 
of superfluid ^He As noted by several authors, the Landau theory is also valid for 
Bose-condensed gases. The thermal conductivity and shear viscosity were first derived for a 
uniform Bose-condensed gas at finite temperatures in a pioneering paper by Kirkpatrick and 
Dorfman ||^. Their results were used by Gay and Griffin to evaluate the temperature- 



dependent damping of first and second sound as predicted by the Landau two-fluid hydro- 
dynamic equations. Our present results are consistent with both of these early papers in the 
appropriate Landau limit, namely when cjr^ ^ 1 However it is important to point out 

that our generalized two-fluid hydrodynamic equations provide a more complete description 
than the original Landau version since they can be used in situations in which the superfluid 
and normal fluid are not in local diffusive equilibrium with each other. 

To illustrate the physics, we use our ZGN' equations in Section IV to study the hydro- 
dynamic normal mode spectrum of a uniform Bose gas in the presence of hydrodynamic 
dissipation. In particular, we show how first and second sound modes are affected by vis- 
cosity and thermal conduction, and also discuss how the new relaxational mode exhibited 
in Refs. is modified. In another paper, we apply these same equations to a discussion 
of the damping of the out-of-phase dipole mode recently observed [l5|,[T2[| in a trapped Bose 
gas. 

II. A REVIEW OF THE ZGN' EQUATIONS 

In this section, we first briefly review the finite temperature ZGN' equations ||^,^ based 
on the assumption that the non-condensate is in local equilibrium. In the next section, we 
calculate the corrections to these equations which arise from a small deviation from local 
equilibrium. The non-condensate atoms are described by the distribution function /(r, p, t), 
which obeys the quantum kinetic equation (we set h = 1 throughout this paper): 

+ P ■ V/(r, p, t)-VU- Vp/(r, p, t) 
at m 

= C,,\f]+C,,[f]. (1) 

Here the effective potential U (r, t) = Ucxt{'r) + 2g[nc{r, t) +n(r, t)] includes the self-consistent 
Hartree-Fock (HF) mean field, and as usual, we treat the inter-atomic interaction in the s- 
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wave approximation with g = Aira/m. The condensate density is ric (r,t) = |$(r,t)p and 
the non-condensate density n(r,t) is given by 

n{r,t) = /(||l/(r,P,t)- (2) 
The two colhsion terms in (|l]) are given by 




x5(p + p2 - P3 - Pi)S{ep + Sp., - Ep.^ - £pj 

X [(1 + + f2)fsU - //2(1 + /3)(1 + /4)] , (3) 



Ci2[f] = 4:TTg^nc J J dp2 J dp3 

x6{mvc + Pi - P2 - P3)S{ec + Sp^ - Sp^ - Ep.^) 
X [5{p - pi) - 5(p - P2) - 5(p - Ps)] 

x[(l + /l)/2/3-/l(l + /2)(l + /3)], 



(4) 



l^r, t) + \mv 



I) takes into account the fact that 
^''r,t) and momentum mVc, 



with / = /(r,p,t), fi = /(r,Pi,t). The expression in 
a condensate atom locally has energy Sciv.t) = fic\^,uj , ^ 
where the condensate chemical potential Uc and velocity will be defined shortly. On the 
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other hand, a non-condensate atom locally has the HF energy ep{r,t) = ^ + U{T',t). This 
particle-like dispersion relation limits our analysis to finite temperatures. 

The equation of motion for the condensate was derived in Ref. |0 (see also Ref. [0) and 
is given by a generalized Gross-Pitaevskii equation for the macroscopic wavefunction $(r,t) 



dt 



-— + f/ext(r) + gn^{r, t) + 2gfi{r, t) - iR{r, t) 
2m 



(5) 



where 



R{v,t) 



with 
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dp 



ri2(r,t) 



■Ci2[/(r,p,t)]. 



(6) 



(7) 



The dissipative term R in (5) is associated with the exchange of atoms between the con- 
densate and non-condensate as described by the collision integral Ci2[/] in (4). We see 
that (1) and (5) must be solved self- consistently. It is customary to rewrite the GP equa- 
tion (5) in terms of the amplitude and phase of $(r, t) = ^JnjT~t)e^^^'^'^\ which leads to 
(v, = V9{r,t)/m) 



dric 
dt 



-ri2[/], 
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m 



dt 

where the condensate chemical potential is given by 



d \ 

+ V, ■ V V, = - V/ic , (8) 



V%Mc(r,t) 

/ic(r, t) = ^ + f/ext(r) + gn,{r, t) + 2gfi{v, t) . (9) 

2mJnc{r,t) 



One sees that in (8) plays the role of a "source function" in the continuity equation 
for the condensate, arising from the fact that C12 collisions do not conserve the number 
of condensate atoms 0. Because of the structure of the equations in (8), they are often 
referred to as "hydrodynamic equations" , even though they are completely equivalent to the 
generalized GP equation in (5). 



Following the standard procedure in the classical kinetic theory of gases [Tj], we take 
moments of (1) to derive the most general form of hydrodynamic equations for the non- 
condensate. These moment equations take the form (// and v are Cartesian components): 

^ + V-(nv„) = ri2[/], 

+ Vn ■ V j Vn^, = - - m{Vnf, - Vc^)Tu[f] , (10) 



di 
di 



+ V ■ (ev„) = - V ■ Q - D^,P^, + 



-m(v„ - Vc)^ + fic-U 



ri2[/]. 



The non-condensate density was defined earlier in (|]) while the non-condensate local velocity 
is defined by 



n(r,t)v„(r,t) = /^^/(r,p,t). (11) 

In addition, we have 





'r, t) = m 


Q 




e 





Vnt,]\ Vnu /(r,p,t). 



^ (p-mv„)2fP - vj/(r,p,t), (12) 



(27r)'^ 2m \m 

(p-mv„)V(r,P,t) 



(27r)3 2m 

Finally, the symmetric rate-of-strain tensor appearing in (10) is defined as 

Formally, these results are exact consequences of the kinetic equation (1). 

The lowest order approximate solution of (1) is based on the assumption that C22 col- 
lisions are sufficiently rapid to force the distribution function to have the form of the local 
equilibrium Bose distribution 
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/(r,P,t) 



(14) 



Here, the temperature parameter f3, normal fluid velocity v„, chemical potential fi, and mean 
field U are all functions of r and t. One may immediately verify that / satisfies C22[f] = 
independent of the value of fl. In contrast, one finds that Ci2[/] is in general finite, namely 

Ci2[/] = 47r^^n,[l - e-^(^-i™(^"-^-)'-^=)] 
dp 



X 



/ (2^ / I ^P3^*^"^^^ + Pi - P2 - P3)S{ei + ec-e2- £3) 



X[6{p- pi) - 6{p - P2) - 5(P - P3)](l + /l)/2/3- 



(15) 



Using the local distribution function (14) to evaluate the moments in (2) and (12), we find 
that the heat current Q(r, t) = 0, and that 



n(r, t) 



dp 



f{r,P,t) 



=0 A3 



93/2[Z) 



(16) 



dp 

(27r)3 3m 



/>,P,t) 



Vn=0 



5, 



(3K- 



(17) 



Here z{r,t) = e^t^-^^'"'*)! is the local fugacity, A(r,t) = [2n /mkBT{r,t)]^/^ is the local 
thermal de Broglie wavelength and gn{z) = Z^i^i-^^V^" the Bose-Einstein functions. The 
kinetic energy density is given by i{r,t) = |P(r, t) which is the same relation as found for 
a uniform ideal gas. 

To summarize, using / ~ /, we obtain the ZGN' lowest-order hydrodynamic equations 
for the non-condensate given in Refs. 



^ + V.(nv„: 



ri2[/], 



-VP - nVU - m(u 



Vc,)Tu[f] 



dP 
'dt 



+ V ■ (Pv, 



2 ~ 2 
■-PV ■ v„ + - 

3 3 



-m(v„ - Vc)^ + ^c-U 



ri2[/]. 



(18) 



where ri2[/] is obtained from (7) with Ci2[/] given by (|15|). 



III. ZGN' EQUATIONS WITH HYDRODYNAMIC DISSIPATION 

We next derive the additional terms which arise from the equations in (10) due to a 
deviation of the distribution function from local equilibrium, / ~ / + 5/ [Q. Following 
Refs. we write this deviation in the form 

5f = /(r, p, t) [1 + /(r, p, tmr, p, t). (19) 

To first order in ip, the C22 collision integral in (3) reduces to 
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C22[/ + 5/]^47r(72 



dp2 f dp 



i2n) 



(2vr)= 



dpA^ip + P2 - P3 - P^jSiei + £2 -£3- £4) 



X//2(1 + /3)(1 + fA){^3 +^4-^2-i^)=L 



22 



(20) 



In the left hand side of (1) and in the C12 colhsion integral, we approximate / by the local 
Bose distribution /. The various derivatives of v„, ft, T and U with respect to r and t 
can be written using the lowest-order hydrodynamic equations given in (18). The resulting 
linearized equation which determines ijj is (for details, see the Appendix) 



u- VT 

f 



mu'^ _ 5g5/2{z) 

2k^T 2(?3/2(^) 
2 



/ mu m 



u ■ w 



ri2[/] 



1 



n 



f(l + f)-Cu[f]=L 



22 



(21) 



Here the thermal velocity u is defined by mu = p— mv„ and w = Vc — v„. The dimensionless 
thermodynamic functions ai, (T2 are defined by 



0-1 (r,t) = 



^mw^ + fic — U 



a2ir,t)=P 



iP7 



U 



|P7 - |n2 



(22) 



where 7(r,t) = ^(7i/2(-2(r, t)). We note that Refs. introduce a related dimensionless 
quantity 7 = (77. 

Since (21) is a linear equation for one may write the solution as = ^/''•^^ where 
if)'^^^ is the solution of 



u- VT 

T 



%5/2(^ 



2/i;BT 25(3/2(2:)] fceT 



m 



-5,,u'\\f(l + f) = L22VlJ 



(1)1 



(23) 



and t/'*-^-* is the solution of 

,2 



mu m 
(T2 H r^'^i + ^ — ' ^ 



3kuT 



^i^/(l + /)-Ci2[/]=L22[^(^)] 



n 



(24) 



We note from its definition that /(r, p, t) is a function of the variable v? while the linearized 
operator L22 defined in (20) is a function of u. The expression (15) for Ci2[/] can be written 
in the form 



dui / d\i2 / (iu3 



m 



Ur, 



X5(w + Ui — U2 — U3)5(/ic — U + -^{w"^ + ""l 
X[5(u-Ui) -(5(U-U2) -(5(u-U3)](l + /i)/2/3, 



(25) 



where /idifr = ji — Hc- This expression shows that Ci2[/] depends on u and w, and in 
particular, obeys the relation Ci2[u, w] = Ci2[— u, — w]. Finally, we note that ri2[/] in (24) 
is independent of u, but depends on w. 
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The equations (23) and (24) can be shown to have unique solutions if we impose the 
following constraints on ip'^^^ and 



dp pjii + m 



dp p^j{i + M'^ = Q. 



(26) 



Physically these constraints mean that the deviations from local equilibrium make no con- 
tribution to fi, v„ and the diagonal component of P^y. The solution for ■j/'*^^^ of (23) has been 
given already in Ref. namely 



(1) 



VT-u 



T 



A{u) + 2D^y u^Uy - -u B{u) 



(27) 



Using this solution for one finds that the heat current density Q and the pressure tensor 
P^y are given by 



S^^uP - 2v 



Q = -kVT + Q 



(2) 



p(2) 



(28) 



where Pj^l^ and Q'-^-* are the contribution from The explicit expressions for the trans- 
port coefficients t] and k are associated with tp^^\ They are given in Eqs. (38) and (34), 
respectively, of Ref. W^ for a trapped Bose gas below Tbec- The analogous transport coef- 
ficients for a uniform degenerate Bose gas above Tbec were first calculated by Uehling and 
Uhlenbeck ||T5 |. 

The additional corrections due to tp^'^^ are given by 



p(2) 
fiu 



m 



dp 



■u^uj{l + f)il) 



(2) 



Q 



(2) 



(27r)3 ^ 



dp m 



(27r)3 2 

The most general solution of (24) for ip^'^^ is of the form 



(29) 



(30) 



If w = v„ — Vc = 0, the left hand side of (24) and -D(u) are isotropic functions of u. Using 
this fact in conjunction with the relation Ci2[u, w] = Ci2[— u, — w] as noted below (25), we 
conclude that D must have the following form 



(31) 



D{u, w) ~ Do{u) + Di{u)w ■ u + D2{u)Wf,Wu yUf,Uu - -u 6f„yj + 0{w 

Using this, the additional terms given by (29) have the following form 

1 



P^^J oc [1 - e-'3('^'^'ff-5™"'')] 
Q(2) oc [1 - e^'^^^'^'ff^s™'"')] 



w. 



(32) 



We note that the constraints in (26) imply that 
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Jdu n7(l + />^'^ = 0, 



(33) 



and thus it follows that the isotropic term Dq{u) in (31) makes no contribution to PjiJ in 
(29). 

In summary, we have obtained the following hydrodynamic equations for the non- 
condensate including the normal fluid transport coefficients 



^ + V-(nv„)=ri2[/], 

(d \ dP dU 

I^IT' \ + ■ ^ ] Vn^^ + ^ h U- 

\ot ax,, ox„. 



d 



dx,. 



2r] 



1 



1^ + V • (ev„) + (V ■ v„)P 



D^.--(TrD)<5, 



-m{yn - Vcfuc - U 



^ fjLV 

dx„ ' 

ri2[/] 



+V ■ (kVT) + 2r] 



1 



- -(TrZ})5 



VQ 



(2) 



(34) 



Since Pj^J and Q'^^^ in (32) are at least of second order in the fluctuations in 5/idifr and 5w 
around static equilibrium, these terms can be neglected when discussing the linearized form 
of these hydrodynamic equations. The equivalent "quantum" hydrodynamic equations for 
the condensate are given in (8). 

In closing this Section, it is useful to summarize the logical structure of our analysis. 
The kinetic equation (1) leads to the exact set of equations involving the variables n, v„, 
£, P^y and Q which are deflned in terms of various moments of the distribution function 
/(r,p,t). In addition, the equations in (10) contain the function in (7) associated 
with collisions between condensate and non-condensate atoms. In Refs. the closed 

set of hydrodynamic equations displayed in (18) were derived by making use of the local 
equilibrium distribution function. In the present paper, we have extended this analysis 
to include a small deviation (19) from local equilibrium, following the Chapman- Enskog 
approach. However, we have only included contributions to the function '?/'(r, p,t) [see (19)] 
associated with the linearized C22 collision integral (Z/22 in (20)), which leads to the linear 
integral equations for in (23) and (24). We have omitted any contribution to '?/'(r, p, t) 
associated with the linearized C12 collision integral. 

Our neglect of the latter contribution can be justifled by noting that C12 in (4) is pro- 
portional to the condensate density nc{r,t). Thus the deviations from local equilibrium 
due to the C12 collision integral are relatively unimportant at temperatures close to Tbec- 
However at lower temperatures, the deviations from local equilibrium due to C12 collisions 
will become increasingly important and one can expect corrections to the values of k, t] and 

obtained in the present paper. Such corrections were evaluated in Ref. |^ for a uniform 
Bose-condensed gas, although these authors made the further assumption that the conden- 
sate and non-condensate were in diffusive equilibrium with each other (i.e., 0). As 
expected, these corrections to the values of k and r] are of order n^/n at temperatures close 
to Tbec (see Eq. (26) of the second reference in 0) and hence can be neglected. 
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IV. NORMAL MODES FOR A UNIFORM BOSE GAS 



In this section, we discuss the normal mode solutions of the linearized hydrodynamic 
equations, as given by 

OSti 

= -V ■ {ricoSvc) - STu, 
^'^^^ (35) 



dt 



and 

d6n 



— — + V ■ (rao5v„) = 6Ti2, 
ot 

_ d5vna 95 P dU ^ dSn d . 
■mno—— = dn- ^gn-o^ h ^— <^ 2r] 



c?(5P 5 ~ 2 ~ 2 

__ = __v . (Po5v„,) + -5v„ ■ VPo + if^co - Uo)ST,2 + -V • (kVST). (36) 

As discussed in Refs. the source term 6ri2 is conveniently expressed in terms of the 

fluctuation of the local chemical potential difference fidis = — fJ-c, 

5Ti2 = (37) 

Tl2 

Here we have introduced an equilibrium relaxation time involving collisions between con- 
densate and non-condensate atoms, 

X(5(pi - P2 - P3)5(/Uc + £pi - £p2 - £p3), (38) 

where it is understood that all quantities now pertain to static thermal equilibrium. 

One can now eliminate ST and Sfidm from the above linearized equations using the 
relations (see (A6) and (A4) of Appendix A) 

6fi = -^5P + 5h + 2g(Sn + Sric), 
no Pono 

6T = ^5P - ^5h. (39) 
Po ^0 

The various coefficients appearing are given by 



— , (720-/^0^0 JZ-^ : 



^ _ Polo ^ _ nl 

2-f070 — 2^0 2-^070 — 2^0 

where /Xco — Uq = —grico in the Thomas- Fermi approximation [0. Using (39) and (40), 
we see that our two-fluid hydrodynamic equations in (35) and (36) reduce to five coupled 
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equations (three for the non-condensate and two for the condensate) in the five variables 5h, 
Suc, Svn, ^Vc, and 6P. Thus one expects a total of five normal modes from the extended 
ZGN' equations. 

We now consider the special case of a uniform gas (t/ext = 0) and look for plane-wave 
solutions ~ g«(k r-ajt) jg Convenient to introduce dimensionless variables 

nc = nc/n, n = n/n, 

P = P/nkBTo, t = r/TBEc, X = gn/kBTBEC, (41) 

where Vd = (5A;b^bec/3?7i)^/^ is the sound velocity of a classical gas at T = Tbec- We also 
introduce dimensionless frequency and wavenumber variables 

Co = ujtq, k = kVciTo, (42) 

where Tq^ = cm(16fcB7BEc/^"^)^^^ is the classical gas collision time ||^, evaluated at T = 
Tbec (but with the quantum collision cross-section for bosons, a = Svra^.) Finally we define 
dimensionless transport coefficients and a dimensionless collision time 

R = K/nv^^TokB, V = rj/nvl^niTQ, = tu/tq. (43) 

In terms of the above dimensionless quantities, we obtain the following closed set of equa- 
tions 

luJOUc = — — — oric + kricoOVc z — —or — —on, (44aj 

3 

iu)5vc = — Xk{6nc + 25n), (44b) 
5 

iu5n = —^6nc + kfiQdVn + 6P + 6n, (44c) 

7-12^ ri2no Ti2no 

6 - 3t - - 4f7 - 

icjSvn = — XkiSuc + 5n) —k5P + —k^5vn, (44d) 

5 brio Suq 

3 3 \ Ti2tno no J ^r^t^ 

3 V ri2tno Po J ^ ^ 

We emphasize that our hydrodynamic equations are restricted by the assumption of local 
equilibrium and thus are only valid if a) -C 1, A; <^ 1. 

It is straightforward to solve the coupled set of equations in (44). Above Tbec, they 
reduce to three equations for three variables, which gives rise to two sound modes {±uk) and 
a thermal diffusion mode. Below Tbec, these equations give two first sound modes {±uik), 
two second sound modes (±^2^;) and one relaxational mode. If the transport coefficients rj 
and K are set to zero, these equations are identical to the ZGN' equations as discussed in 
Refs. [^0. Using (44), we can discuss the effect of r] and k on the first and second sound 
modes, as well as on the relaxational mode. In Fig. 1 of Ref. 0, a graph is given of the first 
and second sound velocities vs. temperature. They are almost identical in both the ZGN 
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limit {uJTfj^ ^ 1) and the Landau limit (cur^ ^ 1). Calculations based on the generalized 
ZGN' equations in (44) show essentially no change in the first and second sound velocities 
when the effects of k and r] are included, for k < 0.5. In Figs. 1-3, we illustrate the effect of 
these transport coefficients on the damping F of the first and second sound modes (where 
Ui = uki—iVi) and on the relaxational mode, all as a function of the dimensionless wavevector 
k. These specific results are for a temperature of T/Tbec = 0.9 and gn/k-BT-^-EC = 0.2, where 
one finds tq/t^ = 0.91, R, = 2.41 and f] = 0.34. We note that in the absence of rj and k, only 
second sound is significantly damped through coupling to the relaxational mode. When we 
include rj and k, Fi is large compared to F2. 

Figs. 1-3 only show the damping in a case where ujt^ ^1. It is interesting to see how 
the relaxational mode depends on the value of cur^j. In Fig. 1 of Ref. @], we have plotted 
the temperature dependence of T12, tq and for gn = 0.1/cb^bec- In Fig- 4, we plot the 
damping of the relaxational mode for k = 0.4, as a function of the temperature. We also 
show the corresponding temperature dependence of t^/tq. As discussed in Ref. [Q, our use 
of a particle-like spectrum with a HF mean-field gives rise to a spurious finite value of the 
condensate density rico at Tbec- As a result, is also finite at Tbec- 

In Appendix B, we derive an analytical expression for the width of the mode centered 
at zero frequency. Working to first order in k, 77, l/ri2 and second order in A = gn/k-QT-Q^Q, 
we obtain 

T^^- + -^^k\ (45) 

As shown in Fig. 4, this approximate expression is in good agreement with the result of a 
direct numerical evaluation of Fr from (44) at temperatures close to Tbec (it deviates more 
and more at lower temperatures because of the non-linear dependence of Fr, on k and I/tu). 
The approximate analytic expression in (45) is useful in that it shows clearly that below T^ec, 
the relaxational mode is strongly coupled to thermal diffusion which arises from the deviation 
from local equilibrium of the non-condensate distribution function. Indeed, one may view 
this damped mode as a renormalized version of the well-known thermal diffusion mode found 
above Tbec BUI- In the classical high-temperature limit, we have a^o = 1 and Pq = nksTQ, 



and (45) reduces in this case to Fr = nk /nCp, where the specific heat Cp = 5/cb/2. This 



is the well-known classical gas result for the thermal diffusion mode |ll0| , ^4 |. 

In the ZGN limit (r^ 00), (45) shows that the damping rate of the relaxational 
mode is mainly due to the finite thermal conductivity (if we set k = 0, we obtain the zero 
frequency mode discussed in Section V of Ref. 0). In the opposite Landau limit of the 
ZGN' equations (r^ 0), the width of this mode is dominated by I/tu, associated with 
the C12 collisions which bring about diffusive equilibrium between the condensate and non- 
condensate components (see Fig. 4). Eq. (45) is in complete agreement with the qualitative 
picture sketched at the end of Section V of Ref. 0. The mode spectrum which the coupled 
ZGN' equations in (44) imply below Tbec is different from the usual Landau two-fluid 
hydrodynamics . In the latter theory, the thermal diffusion mode above Tbec does not 

persist as a relaxational mode below Tbec- Rather, the thermal diffusion mode is interpreted 
as being replaced by two damped second sound modes (±^2^ — iT2) below Tbec- 

Fig. 4 also shows that the width of the relaxational mode increases sharply (over and 
above the spurious jump at Tbec noted above) as the temperature passes from above to 
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below Tbec, as a result of the onset of Cu collisions. This increased width of a mode 
peaked at = will be a useful experimental signature of the new physics implied by the 
extended ZGN' hydrodynamic equations derived in the present paper. Of course, all of the 
above analysis says nothing about the relative weight of the first sound, second sound and 
relaxational modes below Tbec- This requires the evaluation of a specific response function, 
such as the dynamic structure factor [|l^]. Further calculations along these lines are in 
progress. 



V. CONCLUSIONS 

In this paper, we have extended our recent derivation of two-fiuid hydrodynamic 
equations (referred to as the ZGN' equations) to include the effects of a small deviation from 
local equilibrium of the non-condensate atoms. This brings in the usual kind of hydrody- 
namic damping due to the thermal conductivity and shear viscosity of the thermal cloud. In 
calculating these transport coefficients we have only taken into account the deviations from 
local equilibrium due to the C22 collision integral but have neglected the contribution coming 
from the C12 collision integral. This limits the validity of the present calculations to the 
vicinity of the transition temperature. The damping due to hydrodynamic dissipation is in 
addition to that due to the equilibration of the condensate and non-condensate components 
already contained in the original ZGN' equations, as discussed in Refs. @J^. 

For illustration, we presented some numerical results for a uniform Bose gas. In this 
case, we find damped first and second sound modes {±uik — zFi, ±U2k — iT2) and a purely 
relaxational mode u = — iFR. The latter mode is not exhibited by the usual two-fiuid 
hydrodynamic equations which assume that the condensate and non-condensate are in local 
thermodynamic equilibrium. The overall effect of t] and k on the ZGN' predictions in Ref. 
is to introduce additional damping of all three modes. Above Tbec, the relaxational and 
second sound modes (±^2^) merge into a single thermal diffusive mode. The relaxational 
mode below Tbec niay therefore be viewed as the superfluid analogue of the normal thermal 
diffusion mode. 

In another paper [|l^, we shall use these ZGN' equations with transport coefficients to 
discuss the damping of hydrodynamic modes in a trapped Bose gas. In general this is a more 
complex situation, since one must contend with the fact that the local equilibrium solution 
becomes invalid as a starting point in the low density tail of the thermal cloud |TT,TH|. 



However, one can derive general expressions for the damping and renormalization of the 
hydrodynamic modes as given by the ZGN equations. In particular, we shall show that 
the damping of the out-of-phase dipole oscillation of the condensate and non-condensate 
components in a trapped Bose-condensed gas is only weakly affected by the non-condensate 
transport coefficients. As a result, the damping of this mode is entirely due to the collisions 
between the condensate and non-condensate atoms which are responsible for bringing these 
two components into diffusive equilibrium. Further experimental studies of this analogue 
of second sound in superfluid "^He would be of great interest. 
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APPENDIX A: 



We briefly sketch the derivation of the kinetic equation in (21). Using (14) in the left 
hand side of (p, one has 



+ P . V- Vf/(r,t)-Vp 

ot m 



/(r,P,t) 



1 (d_ 
mu / d 



P^ 




m 




' d 




dt 


m 



2A;bT2 \dt m 
VU{r,t) 



■ V V 



u 



/(1 + /)- 



(Al) 



Using the expressions for the density n given by (16) and the pressure P in (17), one finds 



3n 'yk-gT 

an = -—=dT H dz , 

2T z 

dp = —=dT H ^—dz . 

2T z 



(A2) 



where 7 is the variable introduced after (22). One may combine to these equations to obtain 



— — = SdP -dn, 

TP n ' 

^ = -^('-^dn-^-^dP], 
z k^T V 27 2n y ' 

where the thermodynamic functions and (J4 are defined by 

P7 



(A3) 



11? 



IP7 - |n2 



(A4) 



Using the lowest-order hydrodynamic equations given in (18), one finds that the equations 
in (A3) reduce to 



2T 

vn- vr + — airi2[/], 



— = -v„ ■ + a^z — - — . 
at n 



(A5) 
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where the local equilibrium thermodynamic functions ai and o"2 are defined in (22). The 
analogous equation for d\n/dt is given directly by (18). Using these results in (Al), one 
finds that it reduces to 



4- + - • V - VC/ ■ V 

at m 



5P 



{1 ( mu' 

T^' \2kBT ~ 2fik^T 
( mv? m \ 



+ 



m 



ri2[/] 



u 



SkBT 



kuT 



n 



u- (u- V)v- — V-v 

/(! + /), 



(A6) 



where we recall u = p/m — v„. This can be rewritten in the form shown on the left hand 
side of (21). 

Using dz/z = — U)dT/T + djl — 2gdn], one obtains from (A3) 



(7 (J 

dp, = —dP + -^dn + 2g{dn + dn^. 
n pn 



(A7) 



Here a[ and are obtained from ai and (72 in (22) with the replacement of |^«(v„ — Vc)^+//c 
by jl. The relations (A2)-(A4) and (A7) are used in Sec. IV to eliminate (5//difr and 5T from 
the linearized hydrodynamic equations in (35) and (36). 



APPENDIX B: 

In this Appendix, we sketch the derivation of the approximate expression for the relax- 
ation rate Fr given in (45). It is convenient to introduce the five-component vector 

^ Sy^^ §n^^ (Bl) 

The linearized hydrodynamic equations in (44) can then be written in the matrix form 



iujy = Ky, 



(B2) 



where the 5x5 matrix K is given by 



1720 "cO 
"•20 ^^eO _|_ (T40kA:^ 



V 



T12tno _ 

-lAA; 

5 

1720 "cO 

ri2ng 

-lAA; 

5 



no 



g"l(l"eO 

fi2no 
fi2tno Po 

-^k 

5no 

CTl(l"cO 



Ti2no 







nok 

3no 






f\2t 



vo 



3n2f 

-lAA; 

5 

n2^ 

-lAA; 

5 



\ 





ficok 
/ 



(B3) 



Although it is difficult to obtain general analytical solutions of the matrix equation (B2), 
we can obtain approximate analytical expression for Fr by making use of the expression for 
the determinant of this matrix at zero frequency: 



deiK = Tl{nl + Tl){nl + Tl)T^. 



(B4) 
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Working to first order in k, t], l/ri2 and second order in A, a direct evaluation of tlie 
determinant of the matrix in (B3) gives 



no 



O"20 



tPn 



5Po V" tPoCTAO J 

2_ 

-nco<y2Q + 2Pocrio " 2Po) 
5 



(B5) 



If we set K, rj = 0, the first (Qi) and second (^2) sound mode frequencies for small wavevec- 
tors are given by the solution of Eq. (93) of Ref. |0. To second order in A, this gives 



Using (B6) and (B5) in (B4), we obtain 

detK 



To 



Xrin /6 



tPn 



where Tn is given by 



(B6) 



(B7) 



ri2 



no 



\Po + 2gnonco + y^lonlo 



i^o7o - 1^0 

5Po I -^/o- - ^ 

-— + -(2ncocr4o - no) 
2no t 



+ 0{X' 



(B8) 



Neglecting the second small term in the square bracket of (B7), we obtain for Fr the simple 
expression 



TrTo = — + 



2cr4oK7 2 



(B9) 



Numerical calculations show that (B7) is very well approximated by (B9). Using the defi- 
nitions of the various dimensionless quantities in (41)- (43), (B9) can be written in the form 
given in (45). 
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FIGURE CAPTIONS 



FIG.l: Damping of first sound in a uniform gas for T/Tbec = 0.9 and gn = 0.2/cbTbec, as a 
function of the dimensionless wavevector defined in (42). The ZGN' results including 
only K, and only t] are also shown. 

FIG. 2: Damping of second sound in a uniform gas for the same parameters as in Fig. 1. 

FIG. 3: Damping of the relaxational mode vs. wavevector (see Fig. 1). The damping due to 
the thermal conductivity is shown for both the ZGN (r^ oo) and ZGN' theories. 
The effect of the shear viscosity is negligible. 

FIG. 4: Damping of the relaxational mode vs. temperature, for gn = O.I/cb^bec and k = 0.4. 
The ZGN' relaxation time is also shown. All results are normalized to the classical 
gas collision time tq defined below (42). The approximate analytical expression for Fr 
in (45) is compared with a direct numerical solution of the linearized equations in (44). 
The discontinuity at Tbec is spurious, being a result of the mean-field approximation 
used for the thermal excitations |Q. 
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